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UNSTEADY  SEPARATION  AND  UNSTEADY  SHEAR  LAYERS 

(Task  2) 

Principal  Investigator:  Ho,  Chih-Ming 


ABSTRACT 

The  process  of  engulfing  fluids  from  two  streams  into  the  shear  region  is  one  of  the  most 
important  technical  problems  in  fluid  mechanics  research.  A  new  entrainment  mechanism, 
unsteady  azimuthal  deformation  was  identified  to  be  more  effective  than  vortex  merging  in  a 
two-dimensional  flow.  This  mechanism  can  make  the  homogeneous  elliptic  jet  entrain  twice 
the  mass  than  a  cm-'lar  jet.  During  this  contract  period,  the  effect  of  temperature  on  the 
evolution  of  the  three-dimensional  jet  was  studied. 

ACCOMPLISHMENTS 

The  results  are  briefly  summarized  here.  The  detailed  information  can  be  found  in  the 
attached  papers  or  in  the  publication  list. 

Application  of  High  Entrainment  of  Elliptic  Jets 

(Ref.  Ho,  Austin  and  Hertzberg,  Proc.  Asian  Fluid  Mech.  Congress) 

The  high  entrainment  of  a  jet  issuing  from  an  elliptic  nozzle  has  been  applied  in  ramjet 
combustor  and  in  supersonic  jets.  The  data  showed  significant  improvements  in  mixing. 

Stability  Analysis  of  an  Elliptic  Jet 

(Ref.  Koshigoe,  Ho  and  Tubis,  AIAA  paper  No.  87-2733 
Koshigoe,  Tubis  and  Ho,  Phys.  Fluids,  1988) 

The  azimuthal  deformation  of  the  vortices  is  produced  by  the  self-induction  of  the  asymmetric 
distribution  of  vorticity  in  a  structure.  The  vortex  distortion  depends  on  the  stability  process. 
Both  the  amplification  rate  and  the  phase  speed  are  functions  of  the  azimuthal  vorticity 
distribution  and  they  also  vary  with  the  temperature.  When  a  temperature  gradient  exists 
across  the  velocity  shear  with  different  thickness  around  the  nozzle,  the  differences  in  the 
growth  rate  and  the  phase  speed  around  the  nozzle  is  amplified.  Hence,  the  deformation  of 
the  vortices  is  accentuated  by  the  temperature  gradient  and  the  entrainment  is  expected  to 
be  modified. 

The  Hot  Jet  Facility 

We  intend  to  operate  the  jet  in  the  speed  region  where  the  compressibility  effect  is  appreciable 
and  in  the  temperature  range  which  is  much  higher  than  the  room  temperature.  Many 
versions  of  the  modification  have  been  carried  out.  We  were  able  to  increase  the  jet  Mach 
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number  to  about  0.2,  230  ft/sec,  by  installing  three  additional  blowers.  The  maximum 
operating  temperature  is  about  500°  F.  The  temperature  in  the  stagnation  and  at  exits  of 
the  heaters  are  monitored  by  thermistors. 

Flow  Visualization 

There  is  an  appreciable  density  fluctuation  in  the  hot  jet.  We  then  can  visualize  the  en¬ 
trainment  process  in  the  hot  elliptic  jet.  A  shadowgraph  and  a  Schlieren  system  were  set 
up  for  this  purpose.  A  spark-type  light  source  with  micro  second  duration  was  used.  The 
frequency  of  the  light  pulse  was  controlled  by  the  trigger  signal  and  could  be  varied  from  a 
single  pulse  up  to  1  kHz.  We  adjusted  the  pulse  rate  to  match  the  frame  rate  of  the  video 
camera  and  took  a  movie  of  the  hot  jets. 

The  shadowgraph  showed  that  the  initial  boundary  layer  is  laminar  and  vortices  roll  up 
at  a  downstream  location.  The  spreading  in  the  pi  .ne  containing  the  minor  axis  spreads 
much  faster  than  in  the  major  axis  plane.  This  is  consistent  with  the  measurements.  We  did 
not  observe  the  side  jets  caused  by  the  absolute  instability  (Bechert  et  al.,  1988)  in  either 
the  still  picture  or  the  video  movie. 


Temperature  Fluctuations 


The  maximum  temperature  of  the  jet  can  reach  500°F.  To  hold  and  operate  a  hot-wire 
in  a  flow  of  this  high  temperature  is  not  a  trivial  job.  The  time  spent  in  solving  this 
problem  was  more  than  we  expected.  We  need  to  find  a  specially  high  temperature  solder 
to  attach  the  platinum  wire  to  the  steel  prongs.  The  probe  body  also  requires  temperature 
resistance  because  part  of  the  probe  is  inside  the  flow.  The  calibration  procedures  have  been 
established.  The  time  constant  of  the  wire  and  the  digital  frequency  compensation  technique 
were  accomplished.  The  frequency  response  was  extended  from  150  Hz  to  3000  Hz. 

The  pdf  of  temperature  fluctuations  was  measured  along  the  major  and  the  minor  axes  at 
several  streamwise  locations.  When  the  pdf  was  compared  at  corresponding  radial  locations 
in  the  major  and  minor  axis  planes,  an  interesting  feature  was  found;  the  pdf  in  the  minor 
axis  plane  was  being  skewed  toward  the  low  temperature  side.  This  result  reconfirms  the  fact 
that  entrainment  of  the  ambient  cold  near  the  minor  axis  region  is  more  pronounced.  The 
standard  deviation  of  the  pdf  indicates  the  level  of  temperature  fluctuations.  The  standard 
deviations  in  the  minor  axis  plane  are  always  higher  than  those  in  the  major  axis  plane. 

KEYNOTE  TALKS  IN  CONFERENCES 

1.  “Fundamental  Aspects  of  Unsteady  Separation”,  AFOSR  Super  Maneuverability  Ini¬ 
tiative  Meeting,  Boulder,  Colorado,  March  11,  1985. 


2.  “Entrainment  and  Turbulence  Structures”,  22nd  JANNAF  Combustion  Meeting,  Pasadena, 
California,  October  11,  1985. 


□ 
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3.  “Evolution  of  Coherent  Structures  and  Small  Scale  Transition  in  Mixing  Layers”,  Work¬ 
shop  on  Computational  Fluid  Mechanics,  Institute  of  Nonlinear  Science,  University  of 
California  at  Davis,  Davis,  California,  July  18,  1986. 
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4.  “Vortex  Evolution  and  Turbulence  Control”,  24th  JANNAF  Combustion  Meeting, 
Monterey,  California,  October  5-9,  1987. 

5.  “Vortex  Dynamics  of  Unsteady  Flows”,  First  National  Fluid  Dynamics  Congress, 
Cincinnati,  Ohio,  July  23-25,  1988. 
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Vortex  deformation  in  elliptic-core  jets  from  the  perspective  of  linear 
instability  analysis 

Shozo  Koshigoe 
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(Received  10  August  1987;  accepted  17  May  1988) 

An  attempt  is  made  to  identify  the  underlying  mechanisms  for  the  deformation  of  coherent 
structures  that  occurs  in  the  initial  state  of  axis  switching  of  elliptic-core  jets.  The  generalized 
shooting  method  is  applied  to  jets  with  elliptic -core  regions  of  constant  flow.  The  analysis 
shows  that  there  are  three  conditions  on  groups  of  eigenmodes  of  elliptic-core  jets  which  are 
necessary  for  the  deformation.  They  are  ( 1 )  proper  localizations  without  excessive 
overlapping;  ( 2 )  sufficiently  large  phase-speed  differences;  and  ( 3 )  comparable  spatial 
amplification  rates.  The  qualitative  behaviors  of  elliptic -core  jets  in  relation  to  these  three 
conditions  are  studied  with  respect  to  independent  and  joint  variations  of  core  eccentricity, 
azimuthal  distribution  of  momentum  thickness,  and  compressibility. 


I.  INTRODUCTION 

Jets  are  generic  configurations  in  many  engineering  de¬ 
vices.  e.g..  ramjet  combustors,  and  they  therefore  have  at¬ 
tracted  attention  from  numerous  researchers.  Understand¬ 
ing  the  entrainment  mechanism  and  developing  control 
techniques  are  the  key  research  goals  in  jet  studies.  More 
than  a  decade  ago,  effort  was  concentrated  on  investigating 
small  random  eddies,  but  success  was  limited.  Since  the  evo¬ 
lution  of  large  coherent  structures1  was  recognized  to  be  the 
mechanism  that  determined  the  mass  transfer  between 
streams,  progress  has  been  very  much  accelerated.  Many 
active-control2-*  as  well  as  passive-control5  techniques  have 
been  developed  for  manipulating  the  entrainment  process. 
The  basic  control  principle  is  to  enhance  the  unsteady  distor¬ 
tion  of  structures  and  thereby  induce  large  amounts  of  mass 
to  move  across  the  shear  layer.  Passive  control  is  especially 
interesting,  because  no  additional  energy  is  required.  In  the 
case  of  the  elliptic  nozzle  used  by  Ho  and  Gutmark,5  the 
entrainment  can  be  approximately  ten  times  higher  than  that 
for  a  two-dimensional  (2-D)  axisymmetric  or  plane  jet. 
Schadow  and  Gutmark6  tested  this  type  of  passive  control  in 
a  ramjet  combustor  and  found  significant  improvement  in 
the  combustion  efficiency. 

Instability  analysis  offers  a  great  deal  of  insight  into  the 
development  of  free  shear  layers.  Many  properties  of  the 
coherent  structures  in  a>2-D  flow  can  be  predicted  from  in¬ 
stability  calculations.  Even  experiments  in  turbulent  mixing 
layers7  have  demonstrated  that  fine  details  of  the  eigenfunc¬ 
tions  can  be  well  determined  from  the  analysis.  Crighton" 
examined  the  instability  of  an  elliptic  vortex-sheet  jet.  How¬ 
ever.  the  vortex-sheet  approximation  is  only  valid  for  low 
wavenumbirs.  In  the  three-dimensional  (3-D)  elliptic-jet 
case,  new  features,  e.g.,  aspect  ratio  and  distribution  of  the 
momentum  thickness,  are  introduced  into  the  problem  and 
provide  a  new  domain  of  applicability  of  instability  calcula¬ 
tions  to  the  design  of  passive  controls.  There  are  four  classes 


of  instability  modes  for  elliptic-core  jets.6  10  Furthermore,  a 
calculation  scheme  more  general  than  that  used  in  2-D 
cases1 1  is  needed  to  handle  the  nonseparable  boundary  value 
problem. 

One  of  the  most  peculiar  aspects  of  elliptic-core  jets  is 
axis  switching,12  which  is  of  interest  from  a  fundamental 
scientific  point  of  view  and  for  its  potential  applications.  The 
initial  development  of  axis  switching  is  characterized  by  the 
deformation  of  coherent  structures.  This  deformation  can  be 
associated  with  the  self-induction  of  the  asymmetrical  distri¬ 
bution  of  vorticity  in  the  jet,  which  causes  differences  in 
rollup  locations.  The  portions  of  the  elliptic-jet  vortical 
structures  near  the  major  axis  roll  up  slightly  further  down¬ 
stream  from  those  near  the  minor  axis.  Thus  the  overall  pat¬ 
tern  of  rolled  up  structure  is  deformed.  This  deformation  in 
the  elliptic  jet  is  clearly  seen  in  the  forced  elliptic  jets  in  the 
paper5  of  Ho  and  Gutmark  (see  Fig.  1 ).  The  analysis  of  the 
complete  evolution  of  axis  switching  is  beyond  the  scope  of 
the  present  work.  However,  without  the  deformation  of  the 
initial  large-scale-structure  development  as  just  described, 
there  would  be  no  axis  switching. 

In  this  paper,  we  analyze  the  spatial  instability  modes  of 
elliptip-core  jets  in  order  to  demonstrate  the  essential  under¬ 
lying  mechanisms  which  cause  the  initial  vortex  deforma¬ 
tion.  The  generalized  shooting  method, 15  based  on  both  inte¬ 
gral  and  differential  equations,  is  used  for  the  calculations 
since  it  has  been  found  to  be  more  practical  to  apply  than  the 
Helmholtz  integral  equation  formulation.9  Some  of  the  for¬ 
mal  and  numerical  aspects  of  calculations  are  reviewed  in 
Sec.  II  for  convenience  of  discussion  in  the  later  sections. 
(For  more  details,  see  Ref.  13.)  In  Sec.  III.  the  instability 
modes  for  incompressible  flows  are  studied  as  the  limit  ot 
compressible  flows  in  order  to  clearly  assess  the  effects  of 
compressibility  later  in  Sec.  IV. 

Three  conditions  on  groups  of  eigenmodes  are  identified 
as  necessary  for  vortex  deformation  in  elliptic-core  jets. 
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FIG.  1.  Initial  development  region  (»)  and  subsequent  axis-switching  re¬ 
gion  ( s )  of  a  forced  elliptic  jet  in  the  major  axis  (a)  and  minor  axis  (b) 
planes. 


These  conditions' are  studied  in  Sec.  Ill  in  relation  to  the  jet- 
core  eccentricity  and  azimuthal  distribution  of  momentum 
thickness,  and  in  Sec.  IV  in  relation  to  jet  compressibility.  In 
Sec.  V.  streak-line  patterns  are  plotted  for  flow  regimes  rang¬ 
ing  from  low  subsonic  to  supersonic  in  order  to  demonstrate 
explicitly  the  initial  deformation. 


II.  STABILITY  ANALYSIS  FOR  ELLIPTIC  JETS11 

We  express  the  velocity  profile  of  a  3-D  straight  jet  in  an 
inviscid  compressible  fluid  with  incremental  pressure  p  and 
density  p  as 

U  =  W(x,y)et  +  u,  (1) 

where  We.  is  the  jet  mean  flow  velocity  I  W>  0)  and  u  is  the 
fluctuating  velocity  associated  with  spatial  instabilities.  Ap¬ 
plication  of  the  linearized  momentum  and  mass-conserva¬ 
tion  equations  and  the  assumed  (z,r)  dependence, 

(2) 

yields 

(V2  -alc)p  +  Vp-V\n{p'/r/p(a>-aW)z]  =0,  (3) 

where 

a2c  =a:{l  —  [(co/a  —  W)/a]2}.  (4) 

p  is  the  time  averaged  pressure  ( assumed  to  be  independent 
of  position ) , 


V  =  ~e„  — 


.  a 

■  e.. 


dx  +  "ydy' 

and  a  is  the  local  sound  speed,  which  is  related  to  the  time 
averaged  temperature  T  and  density  p  by 


a  =  aoV  T /T0  =  a^pn/p,  ( 6 ) 

where  af)>  Tn,  and  pn  are  the  sound  speed,  temperature,  and 
flow  density  at  the  origin  of  the  coordinates,  respectively.  It 
is  assumed  that  the  temperature  and  the  mean  velocity  are 
related  according  to  the  Busemann-Crocco  law,14 


t=  rjr. 


1 

2 


Ml 


where  T,  =  Tt/T,„  Tt  =  T( v xr  y -  *  ),  Mn  =  W,/a0, 
and  y  is  the  ratio  of  the  specific  heats. 

For  a  given  <s>,  the  solution  of  Eq.  ( 3 )  with  p  finite  and 
—  0  for  -4X-  r  r—  oo ,  yields  the  corresponding  eigenvalue 
wavenumber  ( Re-fa} )  and  spatial  growth  rate  (  —  Im{a} ) . 

We  assume  elliptic-core  jets  whose  mean  flows  W(  r.tb ) 
( r.d>:  polar  coordinates )  have  the  symmetry  properties 

W(r,4  -r  ir)  =  W(r.d)  =  W(r,ir  —  6).  (8) 


The  general  classification  of  eigenmodes  for  this  type  of 
mean  flow  is  discussed  in  Refs.  9-11.  In  this  paper,  we  will 
consider  in  detail  only  the  behavior  of  the  eigenmodes  of  the 
+  +  or  celm  symmetry  class,  for  which 

p(r,ir  —  <b)  —p(r,6)  =  p(r,rr  +  <b).  (9) 

For  circular  jets,  this  class  corresponds  to  the  individual  <b 
dependences  cos  2 mi.  m  —  0, 1 ,2,... .  In  general,  for  jets  with 
mean  flows  satisfying  Eq.  (8),  the  +  +  eigenmodes  may  be 
represented  as 


—  X  :m’(,,)cos  2/naS,  n  =  1,2,3 . 

m  *»  0 

(10) 

Simple  orthogonality  relations  among  the  eigenmodes 
of  Eq.  (3)  do  not  exist.  However,  we  can  apply  the  biortho- 
gonal  relations 15-1 8  derived  in  the  Appendix  to  extract  a  par¬ 
ticular  eigenmode  from  a  general  +  +  expansion  of  p. 

The  spatial  coordinates,  D  and  6,  are  used  to  specify  the 
mean  flow  are  indicated  in  Fig.  2.  The  core  region,  with 
W  =  Wg  constant,  is  bounded  by  the  ellipse  shown.  The  Car¬ 
tesian  coordinates  (x,,yr )  of  a  point  on  this  ellipse  are 


\  Y 

t 


FIG.  2.  Spatial  coordinates  used  to  specify  the  mean  (low  IFf  D.d).  as  de¬ 
scribed  in  the  text.  The  crosshatched  region  in  the  first  quadrant,  where 
\VW\  is  non-negiigible.  is  the  effective  calculation*!  domain  for  Eq.  ( 3 ) . 
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FIG.  3.  Contour  plots  of  the  magnitude  and 
phase  of  the  eigenfunctions  for  the  stream- 
wise-fluctuation-velocity  component  u.,  for 
frequency  a>  =  5.4,  A  / 8  =  2.  and  specifica¬ 
tion  0.02/0.02  for  the  momentum  thickness 
variation  [  the  first  and  second  numbers  refer 
to  the  momentum  thickness  along  the  major 
and  minor  axes,  respectively;  t)„  and  c,  are 
found  from  Eq.  ( 14) ).  The  magnitudes  ( (a) 
and  (c ) }  and  phases  [(b)  and  (d)j  corre¬ 
spond  to  the  most  unstable  eigenmode  of  the 
+  +•  and  -  -  class,  respectively.  The  solid 
and  dashed  lines  represent  positive  and  nega¬ 
tive  phase  contours. 


xt  =  v>0  cos  9, 

-  (11) 
y,  =  ( l/vfo)sin  9, 

where  rn  =  A  /B.  with  A  and  B  being  the  lengths  of  the  semi- 
major  axis  and  semimmor  axis,  respectively;  9  is  the  elliptic 
coordinate  angle;  and  xe  and  y,  are  nondimensionaiized  by 
L  =  mAB  A  given  point  x  =  (x,y)  outside  this  ellipse  is  then 
specified  by  D.  the  distance  from  the  point  to  the  ellipse 
along  a  line  perpendicular  to  the  ellipse,  and  9,  the  elliptic 
coordinate  angle  of  the  point  on  the  ellipse  intersected  by  this 
line. 

By  using  the  symmetry  relations  of  Eq.  (9),  the  compu¬ 
tation  region  for  Eq.  (3)  may  be  effectively  reduced  to  that 
contained  in  the  first  quadrant  (crosshatched  in  Fig.  2).  The 
same  reduction  of  the  computation  region  holds  for  the  other 
symmetry  classes. 1 3 

Calculations  were  done  using  the  mean  flow 

Wtfl,«-±[l-.»h(f±£^)].  (.2) 

where  ©(  9)  is  the  momentum  thickness  and  D,  ( 9)  is  given 
by 


0,(0)  =©(0)  ln((l (13) 

In  Eq.  (13)  fVt  is  a  number  close  to  unity.  All  velocities, 
lengths,  and  temperature  are  nondimensionaiized  by  lVn,  L, 
and  To,  respectively.  For  the  calculations  of  this  paper,  we 
have  used  the  momentum  thickness  function, 

e(0)»0oll+<?,S(0)/S(ir/2>l,  (14) 

where  c,  is  a  constant  for  controlling  the  azimuthal  variation 
of  momentum  thickness  and  S(9)  is  the  arclength  along  the 
curve  c  from  (x,  (0) ,  yt  (0) )  to  (x,  ( 9) ,  y,  ( 9) ). 

In  applying  finite  difference  techniques  to  the  solution 
of  Eq.  ( 3 ) ,  the  crosshatched  region  of  the  x-y  plane  in  Fig.  2 
is  thus  transformed  to  a  rectangular  region  of  the  9-D  plane. 

III.  INSTABILITY  MODES  FOR  INCOMPRESSIBLE 
FLOWS 

In  order  to  assess  the  effects  of  compressibility,  we  first 
investigate  the  behavior  of  the  spatial  instability  eigenmodes 
when  the  flow  is  considered  to  be  incompressible. 

In  Figs.  3  and  4,  contour  plots  are  shown  for  magnitude 
|uz|  and  phase  <bt  of  the  streamwise-velocity-fluctuation 


U-VELOCITY  CONTOURS  U-VELOCITY  CONTOURS 


FIG.  4.  Contour  plots  for  the  magnitude  and 
phase  of  u.  for  the  same  conditions  as  in  Fig. 
3.  The  magnitudes  ((a)  and  (c)  |  and  phases 
((b)  and  (d) )  correspond  to  the  most  unsta¬ 
ble  eigenmodes  of  the  —  +  and  +  —  class, 
respectively.  Most  unstable  eigenmodes  of 
various  classes  are  implied  for  all  subsequent 
figures,  unless  otherwise  stated. 
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FIG.  3.  Contour  plots  of  the  magnitude  and 
phase  of  the  eigenfunctions  for  the  stream- 
wise-fiuctuation- velocity  component  u..  for 
frequency  to  =  5.4.  A  /B  =  3.  and  specifica¬ 
tion  0.02/0.02  for  the  momentum  thickness 
variation  ( the  first  and  second  numoers  refer 
to  the  momentum  thickness  along  the  major 
and  minor  axes,  respectively;  t)„  and  c,  are 
found  from  Eq.  ( 14)  |.  The  magnitudes  [  l  a i 
and  (c )  |  and  phases  [(b)  and  (d)]  corre¬ 
spond  to  the  most  unstable  eigenmode  of  the 
+•  and  —  —  class,  respectively  The  solid 
and  dashed  lines  represent  positive  and  nega¬ 
tive  phase  contours. 


Xt  =  \ro  COS  9, 

_  ■  „  01) 
ye  -  ( 1/V0)sin  9 , 

where  r0  =  A  /B,  with  A  and  B  being  the  lengths  of  the  semi¬ 
major  axis  and  semiminor  axis,  respectively;  9  is  the  elliptic 
coordinate  angle;  and  xe  and  yt  are  nondimensionalized  by 
L  =  mAB  .  A  given  point  x  =  (x,y)  outside  this  ellipse  is  then 
specified  by  D,  the  distance  from  the  point  to  the  ellipse 
along  a  line  perpendicular  to  the  ellipse,  and  9 ,  the  elliptic 
coordinate  angle  of  the  point  on  the  ellipse  intersected  by  this 
line. 

By  using  the  symmetry  relations  of  Eq.  (9),  the  compu¬ 
tation  region  for  Eq.  (3)  may  be  effectively  reduced  to  that 
contained  in  the  first  quadrant  (crosshatched  in  Fig.  2).  The 
same  reduction  of  the  computation  region  holds  for  the  other 
symmetry  classes.13 

Calculations  were  done  using  the  mean  flow 

W(D.Q )  =  —  1  -  tanh  ( -■-—)]  ,  (12) 

2  V  20(0)  /J 

where  ©(  9)  is  the  momentum  thickness  and  D,{9)  is  given 
by 


0J(0)=0(0)ln((l-»')/IFl].  (13) 

In  Eq.  (13)  IF,  is  a  number  close  to  unity.  All  velocities, 
lengths,  and  temperature  are  nondimensionalized  by  IV,,,  L, 
and  To,  respectively.  For  the  calculations  of  this  paper,  we 
have  used  the  momentum  thickness  function. 

(14) 

where  c ,  is  a  constant  for  controlling  the  azimuthal  variation 
of  momentum  thickness  and  S(  9)  is  the  arclength  along  the 
curve  c  from  (xf  ( 0 ) ,  y,  ( 0 ) )  to  (x,  ( 9) ,  y,  ( 9) ). 

In  applying  finite  difference  techniques  to  the  solution 
ofEq.  (3),  the  crosshatched  region  ofthex-y  plane  in  Fig.  2 
is  thus  transformed  to  a  rectangular  region  of  the  9-D  plane. 

III.  INSTABILITY  MOOES  FOR  INCOMPRESSIBLE 
FLOWS 

In  order  to  assess  the  effects  of  compressibility,  we  first 
investigate  the  behavior  of  the  spatial  instability  eigenmodes 
when  the  flow  is  considered  to  be  incompressible. 

In  Figs.  3  and  4,  contour  plots  are  shown  for  magnitude 
|uz|  and  phase  <bz  of  the  streamwise-velocity-fluctuation 
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I 

FIG.  4.  Contour  plots  for  the  magnitude  and 
phase  of  u.  for  the  same  conditions  as  in  Fig. 
3.  The  magnitudes  ((a)  and  (c )  |  and  phases 
[  (b)  and  (d)  1  correspond  to  the  most  unsta¬ 
ble  eigenmodes  of  the  —  +  and  4-  —  class, 
respectively.  Most  unstable  eigenmodes  of 
various  classes  are  implied  for  all  subsequent 
figures,  unless  otherwise  stated. 
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FIG.  6  Core  eccentricity  dependence  of  the  eigenvalues  which  track  from 
those  of  the  m  =  0  (circle).  2  (plus),  and  4  (triangle)  modes  of  a  circular- 
core  let  for  ci  =  5  4  and  0.02/0.02.  in  the  negative-of-the-growth  rate  ( a, ). 
phase  speed  ( C. ).  and  wavelength  (/. )  plane.  The  numbers  in  parentheses 
are  for  A  ,'B 

lar-jet  mode,  an  analogy  to  the  situation  for  the  eigenfunc¬ 
tions.  has  phase  speeds  which  are  intermediate  with  respect 
to  the  phase-speed  extremes  of  the  4-  -j-  0  and  +  +  2 
modes.  The  spatial  growth  rates  for  this  mode  are  compara¬ 
ble  to  those  of  the  other  two  modes.  The  essential  behavior  of 
the  three  eigenmodes  can  be  summarized  by  three  condi¬ 
tions:  ( 1 )  proper  localizations  of  eigenfunctions  without  ex¬ 
cessive  overlapping;  (2)  sufficiently  large  differences  in 
phase  speeds;  and  ( 3 )  comparable  spatial  growth  rates.  Ful¬ 
fillment  of  these  three  conditions  by  the  eigenmodes  may  be 
associated  with  the  deformation  in  the  initial  state  of  the  axis 
switching  of  elliptic-core  jets.  This  point  will  be  demonstrat¬ 
ed  in  Sec.  V. 

B.  Effects  of  azimuthal  distribution  of  momentum 
thickness 

The  Kelvm-Helmholtz  instability  is  associated  with  a 
vorticity  maximum  in  the  velocity  profile.  In  a  2-D  shear 
layer,  the  value  of  the  maximum  vorticity  is  a  constant  in  the 
transverse  direction.  However,  this  is  not  true  in  a  3-D  flow. 
The  azimuthal  variation  of  the  momentum  thickness  will 
certainly  influence  the  characteristics  of  the  instability 
modes  and  can  be  used  as  an  extra  flow  control  property  for 
manipulating  the  instabilty  waves. 

We  now  investigate  the  effect  of  azimuthal  variation  of 
the  momentum  thickness  on  the  eigenmode  behaviors.  For 
A  /B  =  1  the  contour  plots  of  Fig.  7  show  that,  for  all  three 
modes,  the  location  of  maximum  streamwise  velocity  fluctu¬ 
ations  coincides  with  the  location  of  the  minimum  thickness. 
In  contrast  to  the  situation  for  elliptic  jets,  there  is  no  indica¬ 
tion  of  polarization  of  eigenmode  activities  among  the  three 
modes.  Furthermore,  the  associated  variations  of  the  eigen¬ 
values  in  Fig.  8  show  that  the  phase-speed  differences  are 


FIG.  7.  Contour  plots  of  In.:  for  the  evolution  of  the  m  =  0  (a)-(  c  j.  m  =  2 
(d)-(f).  and  m  =  4  (g)-(i)  modes  of  a  circular-core  jet  (A  /B  =  1  j  for 
a  =  5.4  and  variation  of  momentum  thickness  specified  by  0.02/0.02. 0.02/ 
0.022,  and  0.02/0.024  from  top  to  bottom 


greatly  reduced  and  the  higher  modes  have  significant  re¬ 
ductions  in  the  spatial  growth  rates.  Therefore,  the  three 
conditions  are  not  satisfied  for  this  case,  and  the  possibility 
of  simulating  the  ellipuc-jet  behavior  ( the  effects  of  core  ec¬ 
centricity)  using  circular-core  jets  with  nonuniform  mo¬ 
mentum  thickness  distributions  is  remote. 


—  5.6  - 


Cp 


FIG.  8  The  evolution  (with  changes  in  the  momentum  thickness  azimuthal 
dependence )  of  the  m  =  0  ( circle ) .  2  ( tnangle ) .  and  4  ( plus )  modes  of  a 
circular-core  jet  ( A/B=  1)  form  =  5  4.  in  the  negative-of-the-growth  ra>r 
(a, ).  phase-speed  ( Cr  )  plane.  The  numbers  in  parentheses  are  the  values 
for  the  percentage  change  in  momentum  thickness  specification  (  10'7,- 
change  —  0.02/0.022.  etc  ) 
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FIG. 4  Contour  plots  of  u.,  for  the  -  ^  0,  -r  -2.and  -  -  4eigenmodes  (from  top  to  bottomof  the  figure)  for0.02/0.024  (a)-(c),  0.02/0.02  (d)-tD, 
0.024/0.02  lg)-‘  U,  and  m  =  5.4,  A  /B  =  2. 


The  effects  of  azimuthal  variation  of  momentum  thick¬ 
ness  on  eigenfunctions  for  A  /B  —  2  are  given  in  Fig.  9.  The 
contour  plots  for  u.  \  are  displayed  for  momentum  thick¬ 
nesses  (0.02, 0.024)  along  the  (major,  minor)  axes  (a).  It  is 
seen  that  the  spatial  concentrations  of  the  eigenfunctions  are 
shifted  toward  the  major  axis  where  the  momentum  thick¬ 
ness  is  minimal.  For  the  reverse  situation,  with  momentum 
thicknesses  ( 0.024, 0.02 )  along  the  ( major,  minor )  axes  ( c ) , 
it  is  found  that  the  eigenfunction  concentrations  are  shifted 
toward  the  minor  axis  [compare  with  the  uniform  momen¬ 
tum  thickness  case  (0.02,  0.02)  (b)  ] .  Thus  the  effects  of 
momentum  thickness  can  be  used  to  modify  the  eigenfunc¬ 
tion  properties  associated  with  the  core  eccentricity.  The  lo¬ 
calization  characteristics  (condition  1)  among  the  eigen- 
modes  still  exist  for  these  momentum  thickness  variations 
( maximum  variation  20% ) .  However,  a  further  increase  in 
the  variation  can  destroy  the  validity  of  condition  1,  with  all 
of  the  significant  eigenfunction  activities  becoming  concen¬ 
trated  near  the  region  of  minimum  momentum  thickness. 
Note  also  that  the  localizations  of  eigenfunctions  with  mo¬ 
mentum  thicknesses  (0.02/0.024)  show  polarizations  more 
distinct  than  those  for  the  reverse  situation  with  the  momen¬ 
tum  thicknesses  (0.024/0.02).  Therefore,  increasing  mo¬ 
mentum  thickness  along  the  major  axis  while  keeping  the 
momentum  thickness  along  the  minor  axis  constant  appears 
to  lead  to  an  invalidation  of  condition  1 .  A  similar  situation 
is  found  in  the  behaviors  of  corresponding  eigenvalues, 
which  will  be  discussed  in  the  following. 

In  Fig.  10,  we  show  the  variation  of  the  eigenvalues  of 
the  modes  for  A  /B  —  2  when  ( major,  minor)  axis  variations 
of  the  momentum  thickness  are  changed  from  (0.02, 0.024) 
to  ( 0.024. 0.02 )  in  5%  increments  at  three  different  frequen¬ 
cies  ( below,  at,  and  above  the  frequency  of  maximum  ampli¬ 
fication  )  The  effects  of  varying  the  momentum  thickness 
are  clearly  different  for  different  frequencies.  However,  it  is 
seen  that  for  all  frequencies  considered,  the  eigenvalues  for 
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FIG.  10.  Eigenvalues  for  the  +■  +0  (solid  line),  +  2  (dotted  line),  and 

+  +  4  f  dashed  line )  eigenmodes  of  a  jet  with  A  /B  =  2  in  the  negative-cf- 
the-growth  rate  (a,),  phase-speed  (Cp)  plane,  for  to  =  3.6  (a),  5.4  (b).  and 
7.2  (c).  The  variations  of  momentum  thicknesses  are  0.02/0.024  ( <} ),  0.02/ 
0.023  (X),  0.02/0.022 (  + ),  0.02/0.02 1( A),  0.02/0.02(0).  0.021/ 
0.02(V>.  0.022/0.02(8),  0.023/0.02(*),  0  024/0.02 (#). 
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FIG.  1 1.  The  negative-of-che-growth-rate  (a, )  and  phase  speed  (C,)  ver¬ 
sus  frequency  a>  for  eigen  modes  corresponding  to  A  /B  =  2;  -t-  +0. 0.02/ 
0.02  (dashed  line);  -  +•  0.  0.02/0.021  (dotted  line);  +  -r  2,  0.02/0.02 
(solid  line);  and  -r-  1-  2.0.02/0.021  (dash-dot  line). 


the  —  4-0(2)  mode  are  most  sensitive  to  variations  of  mo¬ 
mentum  thickness  in  the  neighborhood  of  either  axis.  At  the 
most  amplified  frequency  (b),  it  is  found  that  increasing  the 
momentum  thickness  along  the  minor  axis  and  keeping  the 
momentum  thickness  along  the  major  axis  constant  in¬ 
creases  the  phase-speed  difference  between  the  4-  4-  0  and 
4-4-2  modes,  whereas  the  reverse  situation  for  the  azi¬ 
muthal  momentum  thickness  along  the  major  axis  constant 
increases  the  phase-speed  difference  between  the  4-4-0 
and  4-4-2  modes,  whereas  the  reverse  situation  for  the  azi¬ 
muthal  momentum  thickness  distribution  decreases  the 
phase-speed  difference.  These  variations  in  azimuthal  mo¬ 
mentum  thickness  distribution  increase  the  differences  in 
growth  rates  for  the  eigenmodes,  however.  Increasing  the 
momentum  thickness  along  the  minor  axis  creates  less 
change  in  the  difference  in  the  spatial  growth  rates  among 
the  eigenmodes  than  the  reverse  case  (e.g.,  comparison 
between  the  cases  (0.02/0.024)  and  (0.024/0.02)  shows 
that  the  difference  in  the  spatial  growth  rates  for  the  first 
case  is  about  half  of  that  for  the  second  case  ] .  Thus  increas¬ 
ing  the  momentum  thickness  along  the  major  axis  appears  to 
invalidate  conditions  2  and  3  for  the  deformation.  In  Fig.  1 1, 
the  eigenvalue  variations  are  given  for  a  continuum  of  co 
values  for  momentum  thicknesses  of  (0.02, 0.02),  and  (0.02, 
0.021)  along  the  (major,  minor)  axes.  It  is  apparent  from 
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FIG.  12.  The  negative-of-ihe- 
growth-rate  (a,)  (a)  and  phase 
speed  (C„)  (b)  versus  frequency  <y 
for  eigenmodes  corresponding  to  A  / 
3=  I;  +  1-0,0.02/0.02;  and  Mach 
number  0  ( solid  line ) .  1  (dotted 
line),  1.5  (dashed  line),  and  1.75 
(dash-dot  line). 


Fig.  1 1  that  selection  of  frequency  (e.g.,  choice  of  5  or  6  for 
a)  does  not  alter  the  above  findings  for  the  eigenvalue  beha¬ 
viors  of  the  two  modes,  4-  4-  0  and  4-  4-2  (e.g.,  the  phase- 
speed  difference  between  two  modes  shown  in  Fig.  6  is  very 
nearly  independent  of  frequency). 

IV.  STABILITY  MODES  IN  COMPRESSIBLE  FLOWS 

In  studying  the  effects  of  compressibility,  the  value  y  for 
the  ratio  of  specific  heats  is  assumed  to  be  1.4,  and  we  only 
consider  the  homogeneous  temperature  case,  with  Tr  =  1  in 
Eq.  (7),  and  (uniform)  momentum  thickness  0.02.  In  Fig. 
12,  the  spatial  growth  rates  (  -  a, )  and  phase  speeds  Cp  are 
plotted  as  functions  of  co  for  the  m  =  0  mode  of  a  circular  jet 
for  Mach  numbers  0,  1.0,  1.5,  and  1.75.  Michalke14  investi¬ 
gated  the  effect  of  Mach  number  on  the  instability  of  circular 
jets.  It  is  not  possible  to  quantitatively  compare  our  results 
with  this  jet  since  the  jet  parameter  (the  ratio  of  jet  radius  to 
the  momentum  thickness)  and  the  Mach  numbers  used  in 
the  two  calculations  are  different.  However,  the  effects  of 
Mach  number  on  the  qualitative  behaviors  of  the  spatial 
growth  rates  are  similar.  The  frequencies  for  the  most  ampli¬ 
fied  modes  for  Mn—  1,  1!5,  and  1.75  are  approximately 
equal  to  5.0,  3.8,  and  2.9,  respectively.  In  Fig.  13,  contour 
plots  of  |ut|  are  given  for  elliptic-core  (A/B-  2)  modes 
4-  4-0,  2,  and  4,  with  M0=  1  (a)-(c),  1.5  (d)-(f),  and 
1-75  (g)-(i)  for  co  =  5.0,  3.8,  and  2.9,  respectively.  It  is  in¬ 
teresting  to  see,that,  with  increasing  eccentricity,  the  eigen¬ 
functions  tracking  from  the  m  =  0(2)  circular-core  modes 
become  more  localized  toward  the  major  (minor)  axis. 
These  trends  are  just  the  opposite  of  those  for  incompressible 
flow  [see  Figs.  5(a)-5(f)  and  Figs.  9(d)  and  9(e)  ].  A  suni- 
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F  .G  '  Cimour  plots  of  u.  for - 0.  -  -  - 4eigenmodes  (from  top  to  bottom  in  the  figure)  for  flow  Mach  number^,,  =  1  (a)-(c),  1.5(d)- 

t  '  i  -:i .  jnd  4  8  =  1.  0.02/0.02.  the  frequency  corresponding  to  maximum  amplification  is  used  for  each  Bow  Mach  number. 


lar  switch  in  the  behavior  of  the  corresponding  eigenvalues  is 
seen  in  Fig.  14  The  switching  in  the  behavior  of  the  eigen- 
modes  occurs  approximately  at  Mach  number  0.45  as  shown 
tn  Fig  15  It  is  clearly  shown  in  Fig.  13  that  the  eigenfunc¬ 
tion  localizations  decrease  for  all  three  modes  as  flow  Mach 
number  increases.  This  tendency  to  invalidate  condition  1  is 
combined  with  one  to  invalidate  condition  3  (due  to  the 
dominant  growth  rate  of  the  4-  •+■  4  mode  as  shown  in  Fig. 
14).  These  result  (see  Sec.  V)  in  a  gradual  decrease  of  rollup 
location  difference  in  streak-line  patterns  as  flow  Mach 
number  increases. 

V.  STREAK-UNE  PATTERNS 

Since  the  eigenfunctions  and  the  corresponding  eigen¬ 
values  are  known,  it  is  possible  to  calculate  the  streak-line 
patterns  issued  from  all  azimuthal  locations  in  order  to  dem¬ 
onstrate  the  deformation  seen  in  Fig.  1 .  However,  it  would 
require  .normous  computing  efforts  to  accomplish  this. 
Therefore,  we  do  the  more  manageable  calculation  of  streak¬ 
line  patterns  emanating  from  the  minor  and  major  axis 
planes.  Another  simplification  in  the  calculation  is  due  to  the 
assumption  that  the  in-phase  forcing  jet  selects  a  linear  com¬ 
bination  of  fluctuation-of-velocity  eigenmodes  of  the  •+■  4- 
class.  The  coefficients  of  the  eigenmodes  may  be  obtained 
from  the  bionhogonal  relations  given  in  the  Appendix.  Mi- 
chaike1*  calculated  the  streak-line  patterns  for  a  plane  shear 
layer  with  a  tanh  velocity  profile,  and  found  the  tendency  of 
the  streak  lines  to  roll  up.  Monkewitz:o  calculated  the  vorti- 
city  contours  at  the  locations  where  the  vortex  merging  oc¬ 
curred  by  taking  subharmonic  resonance  into  consideration. 

In  Figs.  16  and  17.  we  give  the  results  of  streak-line 
calculations  with  time  steps  equal  to  a  quarter  of  the  period, 
based  on  superpositions  of  the  previously  discussed  —  4- 
modesof  an  elliptic  jet.  The  projections  of  the  streak  lines  are 


given  in  the  minor  and  major  axes,  respectively.  The  param¬ 
eter  values  used  were  A  /B  =  2,  co  =  5.4,  and  Af0  =  0.  Using 
Eqs.  ( 1 )  and  ( 2 ) ,  we  followed  the  procedure  of  Michalke 19 
and  set  the  maximum  u-component  magnitude  equal  to 
0.0005  at  z  =  0.  It  is  interesting  to  see  how  the  difference  in 
modal  phase  speeds  is  manifested  in  the  resulting  streak-line 
patterns. 

The  location  of  a  particle  in  the  minor  axis  plane  is  gov¬ 
erned  by  the  equations 

^  =  e-U  +  c.c„  (18a) 

dt 

dx 

—  =  e.-U  4-  c.c.,  (18b) 

dt 

where  c.c.  denotes  complex  conjugate  and  the  right-hand 
sides  are  the  y  and  z  components  of  the  total  velocities.  Simi¬ 
lar  equations  are  obtained  for  the  x  and  z  coordinates  of  the 
particle  position  vector.  The  integration  of  the  particle  tra¬ 
jectory  equations  was  done  using  a  fourth-order  Runge- 
Kutta  method,  with  the  calculational  domain  limited  in  the  z 
direction  by  the  requirement  that  the  maximum  of  the!  dis¬ 
turbance  not  exceed  unity.  For  discussion  of  the  extent  of 
validity  of  the  linear  stability  theory  for  streak-line  calcula¬ 
tions.  see  Refs.  19  and  21. 

It  is  clear  from  Figs.  16  and  17  that  the  difference  in 
phase  speeds  is  reflected  in  the  difference  of  the  rollup  loca¬ 
tions,  with  the  rollup  location  in  the  semimajor  plane  being 
further  downstream  from  that  in  the  semiminor  plane.  The 
rollup  locations  in  the  semiminor  and  major  planes  are  at 
approximately  twice"  the  wavelength  of  the  dominant 
modes  [the  dominant  modes  are  shown  in  Figs.  9(d)  and 
9(e),  respectively].  The  difference  in  the  rollup  locations 
persists  in  the  case  of  .Vf„  =  1 ,  as  can  be  seen  in  Figs.  18(a) 
and  18(b).  The  delays  in  the  rollup  locations  relative  to 


251 1 


ahys  Fluids.  Vol.  31 ,  No  9.  September  1988 


Kosnigoe.  Tubts.  and  Ho 


2511 


-3.5|- 

I 


-3.4k 


0-49  0.50  0.51  0  52  0  53  0.54  0.55  0  56 


Co 


- 1  3i - : - I - : - 

0.49  0  50  0  51  0.52  0.53  0.54  0.55  0.56 

Co 


FIG.  14.  Core  eccentricity  dependence  of  the  eigenvalues  which  track  from 
those  of  the  m  =  0  (circle).  2  (cross),  and  4  (triangle)  modes  of  a  circular- 
core jet  with  0.02/0.02,  for  Mach  number  M„  =  1  (a),  !. 5(b).  1.75(c).  in 
the  negative-of-the-growth  rate  (a.),  phase-speed  (C„l  plane.  The 
numbers  in  parentheses  are  the  values  of  A  /B.  The  frequency  corresponding 
ro  maximum  amplification  is  used  for  each  flow  Mach  number. 


Cp 


FIG.  15.  Core  eccentricity  dependence  of  the  eigenvalues  which  track  from 
those  of  the  m  =  0  (circle)  and  2  (cross)  modes  of  a  circular-core  jet  with 
0.O2/0.02  for  Mach  number  M„  =  0.45.  in  the  negative-of-the-growth  rate 
( a, ) .  phase-speed  ( C„ )  plane.  The  numbers  in  parentheses  are  the  values  of 
A  /B.  The  triangles  indicate  that  the  evolution  tracks  intersect  and  are  thus 
indistinguishable  with  respect  to  their  initial  ( 1 )  points. 

those  in  Figs.  16  and  17  are  a  result  of  the  suppression  of  the 
growth  rates  due  to  the  compressibility  of  the  flow.  The 
streak-line  patterns  for  M  =  1.75  are  shown  in  Figs.  18(c) 
and  18(d).  As  was  discussed  in  Sec.  IV,  the  difficulties  in 
satisfying  conditions  1  and  3,  associated  with  the  gradual 
loss  of  eigenfunction  concentrations  and  of  comparable  am¬ 
plification  rates  among  the  eigenmodes  (due  to  the  domina¬ 
tion  of  4-  +  4  mode  growth  rate),  with  increasing  Mach 
number,  result  in  the  progressive  decrease  of  the  rollup  loca¬ 
tion  difference  as  seen  in  Figs.  16-18. 

VI.  CONCLUSIONS 

An  attempt  has  been  made  to  identify  the  underlying 
mechanisms  for  the  deformation  of  coherent  structures  in 
the  initial  stage  of  axis  switching  of  elliptic-core  jets.  The 
generalized  shooting  method13  was  applied  to  the  spatial  in¬ 
stability  analysis.  The  calculationai  results  showed  that  in¬ 
creasing  core  eccentricity  gives  rise  to  three  distinctive  fea¬ 
tures  of  groups  of  instability  modes:  ( 1 )  proper  localization 
without  excessive  overlapping;  (2)  sufficiently  large  differ¬ 
ences  in  phase  speeds;  and  (3)  comparable  amplification 
rates.  These  features  were  then  identified  as  mechanisms  for 
the  deformation  of  coherent  structures  by  means  of  streak- 
line  calculations.  We  thus  consider  them  as  conditions  on  the 
instability  modes  of  elliptic-core  jets,  which  are  generally 
required  if  deformation  of  coherent  structures  is  to  occur. 

The  qualitative  behaviors  of  elliptic-core  jets  in  relation 
to  these  three  conditions  were  studied  with  respect  to  inde¬ 
pendent  and  joint  variations  of  core  eccentricity,  azimuthal 
distribution  of  momentum  thickness,  and  compressibility. 
The  possibility  of  simulating  elliptic-core  jet  behavior  with 
circular-core  jets  having  nonuniform  momentum  thickness 
was  shown  to  be  remote.  It  was  found  that  increasing  the 
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FIG.  18.  Streak-line  patterns  in  the  minor  axis 
plane  for  M„  =  1  (a).  Mn  =  1.75  ic).  and  in 
the  major  axis  plane  for  M„=  1(b)  and  1.75 
( d ) .  of  an  elliptic -core  jet  with  A /B  =  1. 0.02/ 
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Z 


momentum  thickness  near  the  major  axis  while  keeping  the 
momentum  thickness  near  the  minor  axis  the  same  tended  to 
invalidate  the  above-mentioned  conditions.  The  possibility 
of  realization  of  nonaxis  switching  elliptic-core  jets  is  now 
under  experimental  investigation  using  the  findings  m  this 


paper.  Finally,  it  was  shown  that  the  effect  of  the  compress¬ 
ibility  is  to  deteriorate  the  mechanisms  for  the  vortex  defor¬ 
mation.  Without  the  presence  of  the  shock  structures,  our 
analysis  showed  that  the  deformation  will  not  occur  above 
the  flow  Mach  number  1.75.  Therefore  we  estimate  1.75  as 
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It 

the  upper  limit  of  the  flow  Mach  number  for  the  axis  switch¬ 
ing  with  the  eigenmodes  in  the  ■+■  -+-  symmetry  class. 
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APPENDIX:  BIORTHOGONAL  EXPANSIONS  FOR 
SPATIAL  INSTABILITY  EIGENMODES  OF  EQ.  (S)1*-" 

For  the  nth  eigenmode,  the  eigenvalue  a„ ,  Eq.  ( 3 )  may 
be  written  in  the  equivalent  3x3  matrix  form: 

L(<o)\pn)  -M2)\pn),  (Al) 

where  I p")  is  a  column  vector  whose  components,  pl.pl, 
and  pi,  satisfy  the  boundary  conditions,  p1~ 0  for 
v'x'  +  yr  ;  \t  is  the  local  Mach  number  defined  as 
M  —  W{x)Ja(x)\  and  L(co)  is  the  matrix  operator. 


0 

0 


0 

co{  \  —  M1)  2(1 


L»^V2+  (co'-/a)M- (fV/p)Vp-V  -W  -2coMl 

which  is  independent  of  a„.  It  follows  from  Eqs.  (Al)  and  (A2)  that 
pi  =  apl , 

copl  +  2  WVW-Vpl  =  aWpl, 


( A2) 


(A3) 
( A4) 


and 

(wv1  +  —  \f  -  —  Vp-v]  pi  -  Wpl -IM-pl  =  aW(\  -  At1)  pi  (A5) 

V  a  p  ) 

Elimination  of  pi  and  pi  from  Eqs.  ( A3)-(  A5)  then  yields  Eq.  (3)  for  pi . 

The  adjoint  form  of  Eq.  ( Al )  is 

l  ~(<u)|p'")  =amW(\ -Mz)\pm),  (A6) 

where  pm)  is  the  adjoint  column  vector  whose  components,  p".  pT,  and  p”,  satisfy  the  boundary  conditions.  — 0  for 

v  .x-  -  y  -  30 ;  and  L  ~  ( co )  is  the  adjoint  matrix  operator, 

0  0  V'-W+(o2/a)M  +  V[{lV/p)VpY 


L  ~(co)  = 


0  <u(l  -Mz) 

L^(l-.V/:)  -  V2(l  -\fz)VW 


-  w 

—  2ioMz 


( A7) 


r 


It  is  easily  verified  that  the  lhs  of  the  equation, 
{p""L(ai)'pn)  -  ( p"'L  ~(ca)\pm) 

=  (an-an,)(p"'lV(\-M'-)\pm),  ( A8) 

vanishes,  where  (p"\  is  the  row  vector  transpose  of  Ip"),  and 
rpm'Uco)\p’') 

=  V  [  dx  f  dyp’”L,t{(o)p1,  etc.  (A9) 

I.J  ™  1  J  —  »  V  ”  » 

This  follows  from  the  fact  that  the  area  integral  of  the  lhs  of 
Eq.  (A8)  may  be  converted  into  a  (vanishing)  integral  over 
the  infinite  boundary  of  the  region  of  integration.  Thus,  for 
a,  *crm ,  the  second  factor  on  the  rhs  of  Eq.  ( A8 )  vanishes. 
Using  Eqs.  ( A3)-(A7),  this  result  may  be  expressed  as 

j  dx  j"  dyplUWp?)  =0,  (A10) 

where 

aj-a^W  \a-amfv) 


(1 


Mz) 


a«co\VW\-  \ 

(a>-amW)(a>-a„W)-) 
If  we  now  assume  the  general  expansion 


(All) 


p(x)  =  £  a^pl(x). 


(A12) 


we  may  use  Eq.  (A  10)  to  determine  the  expansion  coeffi¬ 
cients  am  as 


iZ^dxr.  „dyl(Wp?)p(x) 

xdxfZ  ,dyUWp^)pZ(x) 


( A13) 
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ABSTRACT  In  engineering  devices  involving  chemical  reaction, 
it  is  important  to  be  able  to  control  the  mass  transfer  between 
two  streams  of  fluids.  Ho  and  Gutmark  [1J  used  a  small  aspect 
ratio  elliptic  jet  which  can  significantly  enhance  the 
entrainment.  This  efficient  passive  control  technique  is  due  to 
the  unsteady  deformation  of  3-D  vortical  structures  caused  by 
self-induction. 

1. Passive  and  Active  Control  of  Free  Shear  Layers 

The  possibility  of  controlling  the  mass  transfer  in  the 
transverse  direction  of  a  shear  layer  is  important  in  improving 
the  efficiency  of  devices  with  chemical  reaction.  In  a 
two-dimensional  flow,  the  coalescense  of  the  vortical 
structures  in  the  shear  layer  was  identified  to  be  the  main 
entrainment  mechanism  [2].  If  the  vortex  merging  is 
prohibited,  the  growth  of  the  mixing  layer  is  stopped  or  even 
decreased  [3].  On  the  other  hand,  if  the  shear  layer  is 
actively  perturbed  by  the  subharmonics  of  the  most  amplified 
frequency,  multiple  vortices  will  merge  simultaneously  [4]. 
The  growth  of  the  shear  layer  is  increased. 

The  spatial  development  of  the  shear  layer  is  extremely 
sensitive  to  the  initial  perturbations  or  the  boundary 
conditions.  The  alternative  approach  of  controlling  the 
enhancement  is  to  change  the  upstream  boundary  condition  of  the 
shear  layer.  Ho  and  Gutmark  [1,  5]  found  that  a  jet  with  an 
elliptic  nozzle  can  entrain  much  more  fluid  than  that  of  a 
circular  or  a  plane  jet  [Fig.  1].  This  passive  control  method 
is  even  more  advantageous  for  engineering  applications,  because 
no  delicate  forcing  arrangement  is  required.  More  importantly, 
they  identified  a  n^w  entrainment  mechanism;  the  self-induction 
of  an  elliptic  vortex  ring  makes  the  structure  switch  its  axis 
orientation,  the  vortex  element  near  the  minor  axis  moves 
outward  and  makes  a  large  amount  of  ambient  fluid  move  into  the 
jet  near  the  minor  axis  region.  This  concept  can  be 
generalized  and  used  in  other  flow  configurations,  such  as 
combustion  chamber. 


2 . Applications  in  Combustion 

The  ramjet  is  a  device  which  can  be  benefited  by  the 
entrianment  contol  technique.  There  is  a  short  distance  from 
the  flame  holder  to  the  exhaust  nozzle.  Combustion  needs  to  be 
accomplished  during  a  short  residence  time  and  combustion 


instability  has  been  a  troublesome  problem.  The  advantages  of 
using  an  elliptic  .jet  in  the  ram.jet  has  been  shown  to  be 
phenomenal.  Schadow  et  al .  [6]  used  a  3:1  elliptical  nozzle 

in  a  jet  with  combustion.  They  found  that  the  centerline 
temperature  increased  sharply  a  short  distance  from  the 
elliptic  nozzle.  At  the  end  of  the  jet  potential  core,  the 
temperature 

was  much  higher  than  that  of  a  circular  jet  [Fig.  2]. 

The  combustion  instability  problem  was  alivated  by  using 
another  type  of  asymmetric  nozzle.  They  used  a  triangular 
nozzle  and  injected  fuel  near  the  tips  of  the  triangle  [7].  In 
this  way,  the  fuel  is  mixed  by  the  small  eddies  near  the  tips 
and  the  large  structures  do  not  trigger  the  combustion 
instability. 

3. Supersonic  Asymmetric  Jet 

The  spreading  rate  of  a  supersonic  shear  layer  is  much 
slower  than  that  of  a  subsonic  flow  and  the  combustion 
efficiency  of  supersonic  flame  is  hindered  by  the  low 
mixing  rate.  Actually,  this  is  the  most  pressing  problem  in 
developing  a  hypersonic  aircraft.  In  a  preliminary  experiment, 
we  found  that  the  small  aspect  ratio  rectangular  nozzle  could 
enhance  the  entrainment  as  it  did  in  the  subsonic  flow  [8]. 

In  general,  the  supersonic  flows  spread  slower  than  that 
of  subsonic  flows.  However,  the  entrainment  improvement  of  the 
rectangular  nozzle  over  the  circular  nozzle  is  about  the  same 
for  both  supersonic  and  subsonic  jets  in  the  far  downstream 
locations.  Near  the  Nozzle,  the  variations  of  the 
cross-section  area  ratios,  indicating  the  shear  layer 
spreading,  are  similar  for  supersonic  and  subsonic  cases. 
Hence,  we  expect  that  the  supersonic  rectangular  jet  will 
entrain  more  mass  than  that  of  a  supersonic  circular  jet  in  the 
region  near  the  nozzle. 

4. Asymmetric  Jet  in  Confinement 

Most  of  the  combustion  is  taking  place  inside  \  confined 
space.  The  entrainment  process  is  very  different  from  a  shear 
layer  in  the  free  space.  As  has  been  pointed  out  in  our 
previous  findings  [1]  the  unsteady  evolution  of  the  vortex 
structures  can  engulf  the  fluids  into  the  mixing  region.  We 
used  a  2:1  aspect  ratio  rectangular  jet  in  a  confined 
environment.  The  deformation  of  the  vortical  structures  was 
found  to  be  much  more  convoluted  than  those  in  a  free  jet. 
These  deformations  are  produced  by  the  local  3-D  shear  layer 
and  the  induction  of  the  image  vortex.  In  other  words,  a 
properly  designed  confinement  should  be  able  to  facilitate  the 
mixing  process. 
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CENTERLINE  TEMPERATURE  OISTWBUTION 


Fig,  1  Mass  entrainment  rate 


Fig.  2  Centerline  temperature 


SEARCH  FOR  CHAOS  IN  FREE  SHEAR  FLOWS 
THE  CASE  OF  THE  WAKE-SHEAR  LAYER 

(Task  3) 

Co-Principal  Investigators:  L.  G.  Redekopp,  P.  Huerre,  and  C.-M.  Ho 


STATUS  OF  RESEARCH  EFFORT 
1  Local  and  Global  Instabilities 

Many  flows  of  technological  interest  such  as  jets,  wakes,  mixing  layers,  boundary  layers, 
etc.  are  spatially- developing  (i.e.,  the  mean  velocity  profile  is  non-uniform  in  the  streamwise 
direction).  In  the  past  three  years,  we  have  been  able  to  develop  a  general  theoretical 
framework  to  describe  these  flows  from  a  hydrodynamic  stability  perspective.  The  status  of 
recent  theoretical  and  experimental  efforts  by  our  group  and  other  teams  has  been  reviewed  in 
[15]  (numbers  in  brackets  refer  to  respective  publications  listed  at  the  end  of  the  report).  Our 
approach  differs  markedly  from  previous  work:  instead  of  examining  the  spatial  development 
of  downstream-going  waves  at  a  given  external  real  frequency,  we  argue  that  one  must  first 
determine  the  absolute/ convective  nature  of  the  local  instability  at  each  streamwise  station. 
If  the  entire  flow  is  locally  convectively  unstable,  as  in  spatial  mixing  layers,  boundary 
layers,  constant-density  jets,  flat-plate  wakes,  etc.,  the  classical  method  mentioned  above  is 
appropriate.  However,  when  a  sufficiently  large  region  of  absolute  instability  exists  within 
the  flow  domain,  as  in  bluff  body  wakes,  low  density  jets,  capillary  jets,  etc.,  one  must 
resort  to  a  global  mode  description  which  includes  upstream  and  downstream-going  waves  of 
unknown  complex  frequency.  The  global  mode  frequency  and  its  spatial  structure  are  then 
obtained  by  solving  an  eigenvalue  problem  in  the  streamwise  direction. 

Simplified  amplitude  evolution  models,  such  as  the  Ginzburg-Landau  equation,  have 
proved  to  be  extremely  useful  to  test  these  new  theoretical  concepts.  We  have  shown  on 
these  models  that  the  onset  of  global  instability  results  from  a  well-defined  sequence  of 
bifurcations  [2,6].  Detailed  stability  analyses  of  wakes  behind  bluff  bodies  have  confirmed 
that  the  same  scenario  is  observed  experimentally  in  bluff-body  wakes  [7,8].  The  effects  of 
forcing  on  global  mode  dynamics  have  also  been  studied  on  the  basis  of  simplified  models: 
it  has  been  demonstrated  that  the  optimal  station  for  the  excitation  of  global  modes  is 
located  well  upstream  of  the  point  of  maximum  absolute  growth  rate  [4,11].  Furthermore, 
the  application  of  a  forced  evolution  model  to  the  description  of  the  preferred  mode  in 
homogeneous  jets  has  led  to  surprisingly  good  results:  the  predicted  Strouhal  number  of  0.225 
compared  very  favorably  with  the  measured  value  of  0.25  [11,3].  This  led  us  to  conclude  that 
the  preferred  mode  in  jets  is  a  weakly-damped  global  mode  which  is  selectively  destabilized 
by  extremely  low  levels  of  external  noise.  Finally,  a  study  of  the  secondary  instabilities 
arising  in  some  amplitude  equations  has  indicated  that  primary  and  secondary  instabilities 
may  have  a  distinct  absolute/convective  character  within  the  same  flow  [1]. 

As  briefly  alluded  to  earlier,  hydrodynamic  stability  calculations  of  wake  profiles  [7,8] 
have  also  been  carried  out  to  determine  whether  absolute/convective  instability  concepts  are 
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pertinent  in  real  flows.  These  studies  have  revealed  that  a  region  of  local  absolute  instability 
exists  in  the  near  wake  behind  bluff  bodies  when  the  Karman  vortex  street  is  self-sustained. 
As  the  Reynolds  number  is  decreased  below  critical,  the  absolutely  unstable  region  persists, 
but  it  is  too  small  to  lead  to  a  self-sustained  global  mode.  In  this  case,  the  flow  may  be  locally 
unstable,  but  there  is  no  spatially-coherent,  temporally-periodic  structure  characteristic  of 
the  Karman  vortex  street. 

A  much  closer  link  between  local  and  global  properties  has  been  established  within  the 
WKBJ  approximation  for  a  general  class  of  partial  differential  equations  in  one  spatial  di¬ 
mension  [9,10].  Rigorous  mathematical  results  have  been  derived.  In  particular,  a  simple 
frequency  selection  criterion  has  been  proposed:  the  global  mode  frequency  is  given  by  the 
value  of  the  absolute  frequency  o>o(A,)  at  the  stationary  point  X,  such  that  ( dw0/dX)x ,  =  0. 
This  provides  an  explicit  relationship  between  global  and  local  characteristics.  It  is  striking 
that  the  global  frequency  only  involves  a  single  streamwise  location  in  the  complex  X-plane. 

2  Spatial  chaos:  model  studies 

In  convectively  unstable  flows,  interactions  between  different  modes  lead  to  a  radically  dis¬ 
tinct  type  of  dynamics.  In  contrast  with  absolutely  unstable  flows,  one  expects  a  great 
sensitivity  to  external  perturbations.  We  have  conducted  an  asymptotic  study  of  the  inter¬ 
action  between  the  sinuous  and  varicose  modes  in  wake-shear  layers  when  the  critical  layer 
is  dominated  by  viscous  effects  [5].  The  complex  amplitudes  of  each  mode  have  been  shown 
to  be  governed  by  a  system  of  coupled  subharmonic  resonance  equations.  A  numerical  in- 
vestigaton  of  the  resulting  dynamical  system  has  revealed  a  great  variety  of  possible  spatial 
vorticity  distributions  as  the  forcing  frequency  is  gradually  increased.  In  all  cases  the  flow 
is  time  periodic  but  spatial  chaos  can  be  generated  for  specific  ranges  of  forcing  frequencies. 
The  nature  of  the  spatial  disorder  is  highly  intermittent:  the  sinuous  mode  dominates  the 
evolution  but  localized  bursts  of  the  varicose  mode  occur  randomly  in  space.  This  inter- 
mittency  has  been  used  to  advantage  to  derive  an  analytical  one-dimensional  map  relating 
successive  peaks  of  the  varicose  mode. 

3  Spatial  chaos:  experimental  studies 
3.1  The  water  channel 

The  wake-shear  layer  water  channel  was  assembled  and  preliminary  testing  occurred  during 
the  current  contract  period.  The  facility  has  a  splitter  plate  which  divides  the  stagnation 
and  contraction  sections  allowing  control  of  the  two  independent  streams  which  merge  down¬ 
stream  of  the  splitter  plate  in  the  test  section.  A  major  modification  of  the  trailing  edge 
of  the  splitter  plate  was  necessary  after  early  tests  revealed  that  local  separation  caused  a 
high  level  of  random  disturbances.  The  original  trailing  edge  was  made  of  fiberglass  and 
the  tapper  angle  was  larger  than  prescribed  in  the  original  design.  We  have  replaced  the 
last  section  of  the  splitter  plate  with  a  trailing  edge  made  of  precision  machined  brass.  The 
local  separation  problem  was  found  to  be  eliminated.  The  speeds  in  the  streams  on  both 
sides  of  the  splitter  plate  can  be  independently  controlled  for  flow  speeds  up  to  1  ft/sec. 
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Suction  or  blowing  from  the  trailing  edge  is  also  possible  to  provide  independent  control  of 
the  wake  component  of  the  velocity  profile  for  the  purpose  of  investigating  various  types  of 
instabilities  and  mode  interactions. 

3.2  Instrumentation 

A  powerful  new  type  of  instrument,  a  Particle  Displacement  Image  Velocimeter  (PDIV), 
will  be  used  to  measure  the  velocity  field.  The  PDIV  instrument  provides  a  measurement 
of  several  thousands  of  velocity  vectors  at  one  time.  With  this  capability,  it  is  possible 
to  investigate  the  spatial  chaos  problem  in  that  the  instantaneous  spatial  structure  of  the 
velocity  field  can  be  obtained  as  a  function  of  time. 

All  components  of  the  PDIV  system  have  been  received.  We  experienced  a  problem 
with  the  SUN  computer  which  was  subsequently  found  to  arise  from  an  incompatible  disk 
controller.  Fortunately,  we  managed  to  resolve  this  problem  in  the  first  week  of  December. 
The  whole  system  is  now  functioning  and  should  be  in  operating  condition  by  January. 

3.3  Experimental  Results 

When  both  streams  have  the  same  speed,  a  wake  flow  is  produced  downstream  of  the  trailing 
edge  of  the  splitter  plate.  The  arrangement  of  coherent  structures  was  observed  to  be  of  the 
Karman  vortex  street  type.  They  were  counter-rotating  and  asymmetrically  placed.  When 
the  velocity  difference  between  the  two  streams  was  increased,  a  mixing  layer  eventually 
formed.  All  the  structures  rotated  in  the  same  direction.  There  must  be  a  velocity-ratio 
range  such  that  both  the  Karman  mode  and  the  mixing-layer  mode  compete.  We  placed  a 
pulsed  hydrogen  bubble  wire  in  the  flow  and  visualized  the  vortex  patterns  from  the  streak 
lines.  When  the  velocity  ratio  is  close  to  zero  (i.e.,  strong  wake  influence),  the  vortices  shifted 
from  the  wake  type  to  the  mixing  layer  type  intermittently.  Whenever,  the  wake  type  vor¬ 
tices  prevailed,  the  thickness  of  the  sheared  region  increased.  A  hot-film  probe  was  used 
to  measure  the  passage  frequency  of  the  vortices.  Based  on  the  time-averaged  spectrum, 
only  one  peak  could  be  identified.  In  other  words,  the  flow  was  temporally  periodic,  but 
the  vortex  arrangement  exhibited  significant  spatial  variability.  We  further  found  that  the 
phenomenon  described  above  only  occurred  over  a  very  narrow  range  of  velocity  ratios. 
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P.  Huerre,  Invited  Lecture,  March  Meeting  of  the  American  Physical  Society,  St.  Louis, 
Missouri,  March  21,  1989. 
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14.  “Nature  of  Bluff-Body  Near- Wake  Instabilities”,  P.  Huerre,  Lead-off  Talk,  ONR  Ac¬ 
celerated  Research  Initiative  on  bluff  body  wake  dynamics,  Lehigh  University,  Penn¬ 
sylvania,  May  23,  1989. 

15.  Open  Flow  Instabilities”,  P.  Huerre,  Invited  Lecture,  Sixth  Taylor- Vortex  Flow  Work¬ 
ing  Party,  Universite  Libre  de  Bruxelles,  Bruxelles,  Belgium,  May  30,  1989. 

PROFESSIONAL  PERSONNEL 

Dr.  C.  M.  Ho:  Co-Principal  Investigator. 

Dr.  P.  Huerre:  Co-Principal  Investigator. 

Dr.  L.  G.  Redekopp:  Co-Principal  Investigator. 

Dr.  M.  Rossi:  Research  Associate. 

Ms.  D.  Wallace:  Research  Assistant. 

Mr.  G.  Vance:  Research  Assistant. 
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